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Abstract 

Recent observations suggest that 7-ray bursts (GRBs) and their afterglows 
are produced by jets of highly relativistic cannonballs (CBs), emitted in su- 
pernova (SN) explosions. The CBs, reheated by their collision with the shell, 
emit radiation that is collimated along their direction of motion and Doppler- 
boosted to the typical few-hundred keV energy of the GRB. Accompanying 
the GRB, there should be an intense burst of neutrinos of a few hundreds of 
GeV energy, made by the decay of charged pious produced in the collisions 
of the CBs with the SN shell . The neutrino beam carries almost all of the 
emitted energy, but is much narrower than the GRB beam and should only 
be detected in coincidence with the small fraction of GRBs whose CBs are 
moving very close to the line of sight. The neutral pious made in the transpar- 
ent outskirts of the SN shell decay into energetic 7-rays (EGRs) of energy of 
0(100) GeV. The EGR beam, whose energy fluence is comparable to that of 
the companion GRB, is as wide as the GRB beam and should be observable, 
in coincidence with GRBs, with existing or planned detectors. We derive in 
detail these predictions of the CB model. 
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1 Introduction 



For a third of a century, gamma ray bursts (GRBs) have constituted a great as- 
trophysical mystery. Their origin is still an unresolved enigma, in spite of recent 
remarkable observations in the field: the discovery of GRB afterglows 0], the 
discovery |^ of the association of GRBs with supernovae (SNe), and the measure- 
ments of the redshifts [Q of their host galaxies. The current generally accepted 
view is that GRBs are generated by synchrotron emission from relativistically ex- 
panding fireballs, or firecones, produced by collapses or mergers of compact stars 
P, by failed supernovae or coUapsars 0, or by hypernova explosions 0. It was 
further suggested that these highly relativistic fireballs produce large fiuxes of very 
high energy neutrinos in coincidence with the GRBs 0. But various observations 
suggest that most GRBs are produced by highly collimated superluminal jets and 
not by relativistically expanding fireballs 0, [1^, [Tl|, |12|. 



In a recent series of papers |jT2|, |22|, |23| we have outlined a cannonball (CB) model 
of GRBs which, we contend, is capable of describing the GRB phenomenology, and 
results in interesting predictions. The CB model is based on the following analogies, 
hypothesis and explicit calculations: 



Jets in astrophysics. Astrophysical systems, such as quasars and microquasars, in 
which periods of intense accretion into a massive object occur, emit highly collimated 

1 / y'l — v^/c^ of these jets ranges from mildly 
to quite relativistic: 7 = C (10) for 



jets of plasma. The Lorentz factor 7 
relativistic: 7 ~ 2.55 for PSR 1915+13 
typical quasars 



H 
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and even to highly relativistic: 7 ~ 10'^ for PKS 0405—385 
These jets are not continuous streams of matter, but consist of individual 
blobs, or "cannonballs" . The mechanism producing these surprisingly energetic and 
collimated emissions is not understood, but it seems to operate pervasively in nature. 
We assume the CBs to be composed of ordinary "baryonic" matter (as opposed to 
e"'"e~ pairs), as is the case in the microaquasar SS 433, from which Ly^, and metal 



Kq lines have been detected 27, 28 



The GRB/SN association. The original observation of a spatial and a temporal 
coincidence between GRB 980425 and the relatively close-by supernova SN 1998bw 
(redshift z = 0.0085), that suggested a physical association [|], has developed into 
a much more convincing case for the claim [0 that many, perhaps all, of the long- 
duration GRBs are associated with SNe. Indeed, of the dozen and a half GRBs 
whose redshift is known, the nearest six that have redshifts z < 1 show in their 
afterglow an additive "bump", with the time dependence and spectrum of a SN 



akin to 1998bw properly corrected 
galactic extinction |f3|: GRB 990228 0, p 
0; GRB 990712 



TBI for the different redshift values and 



T6|; GRB 970508 [Q; GRB 980703 



23; GRB 991208 0]; GRB 000418 |^|. 
In all the other cases with larger redshifts there is one or more good reasons for 
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such a bump not to have been seen: no observations of the afterglow at late time are 
available, the expected bump is below the sensitivity of the late time observations, 
the spectrum of SN 1998bw is not known at the frequencies required to extrapolate 
its light curve to a much higher z. Thus, observationally, seven out of sixteen — and 
perhaps all — of the GRBs of known redshift have a SN associated with them. The 
energy supply in a SN event similar to SN 1998bw is too small to accommodate the 
fluence of cosmological GRBs, unless their 7-rays are highly beamed. SN 1998bw is 
a peculiar supernova, but that may be due to its being observed close to the axis 
of its GRB emission. It is not out of the question that a good fraction -perhaps 
all- of the core-collapse SNe be associated with GRBs. To make the total cosmic 
rate of GRBs and SN compatible, this nearly one-to-one GRB/SN association would 
require beaming into a solid angle that is a fraction f ~ 2 x 10^^ of An [[T^. The 
CB model, for the emission from CBs moving with 7 ~ 10^, implies precisely that 
beaming factor (the numerical details are reproduced in Section 9). 

The GRB engine. We assume a core-collapse SN event not to result only in 
the formation of a central compact object and the expulsion of a supernova shell 
(SNS). A fraction of the parent star's material, external to the newly-born compact 



object, should fall back in a time very roughly of the order of one day |]12[ and, 
given the considerable specific angular momentum of stars, it should settle into an 
accretion disk and/or torus around the compact objectQ The subsequent sudden 
episodes of accretion — occurring with a time sequence that we cannot predict — 
result in the emission of CBs. These emissions last till the reservoir of accreting 
matter is exhausted. The emitted CBs initially expand in the SN rest-system at a 
speed /9c/7, with f3c presumably of the same order or smaller than the speed of 
sound in a relativistic plasma (/? = l/\/3). The solid angle a CB subtends is so 
small that presumably successive CBs do not hit the same point of the outgoing 
SNS, as they catch up with it. These considerations are illustrated in Fig.(|l|). 

The GRB. From this point onwards, the CB model is not based on analogies 
or assumptions, but on processes whose outcome can be approximately worked out 
in an explicit manner. The violent collision of the CB with the SNS heats the CB 
(which is not transparent at this point to 7's from tc^ decays) to a temperature 
that, by the time the CB reaches the transparent outskirts of the SNS, is ~ 150 



eV, further decreasing as the CB travels |2^. The resulting CB surface radiation. 



Doppler-shifted in energy and forward-collimated by the CB's fast motion, gives 
rise to an individual pulse in a GRB, as illustrated in Fig.(|I|). The GRB light curve 
is an ensemble of such pulses, often overlapping one another. The energies of the 
individual GRB 7-rays, as well as their typical total fluences, indicate CB Lorentz 
factors of O(IO^), as the SN/GRB association does p2[. The GRB properties most 



relevant to the current investigation are reviewed in Section 5. 



^We choose to base our conjectures on analogies with known processes, as opposed to computer 
simulations. The latter do not yet realistically include rotation, magnetic fields, the transport of 
angular momentum... Most noticeably, they do not produce SN explosions. 
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The GRB's afterglow. The CBs, after they exit the SNS, cool down by bremsstrahlung 
and radiate by this process, by inverse Compton scattering, and by synchrotron ra- 
diation of their electrons on their enclosed magnetic field, much as the plasmoids 



emitted by quasars and microquasars do |12|. The CB model provides an excellent 



detailed description of optical afterglows [|13[. The early afterglow spectrum and 
light curve are complicated by the fact that, about a day after the GRB emission, 
CBs cool down to a temperature at which e-p recombination into H takes place. 
This gives rise to Ly-a lines that the CB's motion Doppler-shifts to (cosmologically 
redshifted) energies of order a few keV, an energy domain that, interestingly, coin- 
cides with that of the Fe lines that an object at rest would emit. Recombination 
also gives rise to a multiband-flare in the afterglow. These CB-model's expectations 
are in good agreement with incipient data on X-ray lines and flares 

In all of the above considerations we have exploited the GRB/SN association 
to conclude that GRBs, at least the long duration ones, are associated with core- 
collapse SNe. This allows us to be very specific in our predictions []I2|, |23 



concerning the collision of the CBs with the SNS, and the consequent properties of 
the GRB pulses (the density profile of SNSs is known from observations; the typical 
energy of CBs we can infer from the assumption that the large peculiar velocities of 
neutron stars are due to an imbalance between the momenta of the jets of CBs they 



emit as they are born ||TT|). But the sites of GRB emission may not be only SNe. 
Any process of violent accretion, such as a merger between neutron stars or other 
compact objects, may result in the emission of CBs. If the latter encounter matter 
on their way, such as circumstellar gas or a molecular cloud, the processes leading 
to 7-ray emission would be similar to the ones pertaining to a SN engine. 

In this paper we address two other concrete predictions of the Cannonball Model: 
the emission of neutrinos and of energetic 7-rays (EGRs). Once again, to be specific, 
we exploit our explicit model of GRBs emitted in core-collapse SN events. The 
neutrinos are made by the chain decays of charged pions, produced in the collisions 
of the CBs' baryons with those of the SNS, as in Fig.(|l]). The u beam carries almost 
all of the emitted energy, but is much narrower than the GRB beam and should only 
be detected in coincidence with the small fraction of GRBs whose CBs are moving 
extremely close to the line of sight. The EGRs are made by the decay of neutral 
pions, but only from 7r° production close enough to the outskirts of the SNS for 
the 7-rays not to be subsequently absorbed, see Fig.(|I]). The EGR beam, whose 
fluence is comparable to that of the GRB, is as wide as the GRB beam and should 
be observable, in coincidence with GRBs, with existing or planned detectors. The 
EGR beam peaks at energies of tens of GeVs, while the u beam is about one order 
of magnitude more energetic. 
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2 Times and energies 



Let 7 = — f3'^ = Ecb/(Mcbc^) be the Lorentz factor of a CB, which diminishes 

with time as the CB hits the SNS and as it subsequently plows through the inter- 
stellar medium. Four clocks ticking at different paces are relevant to a CB's history. 
Let tsN be the local time in the SN rest system, tcB the time in the CB's rest system, 
tob the time measured by a nearby observer viewing the CB at an angle 6 away from 
its direction of motion, and t the time measured by an earthly observer viewing the 
CB at the same angle, but from a "cosmological" distance (redshift z 7^ 0). Let x 
be the distance traveled by the CB in the SN rest system. The relations between 
the above quantities are: 

dx 

dtsN = 7 dtcB = -5- ; 

p c 

dtcB = S dtob ; 

dt = (1 + z) dtob = dtsN , (1) 

7 



where the Doppler factor 6 is: 



S^^^ ^-1— (2) 

7(l-/?cos^) ~ (l + ^V) ' 



and its approximate expression is valid for ^ 1 and 7^1, the domain of interest 
here. Notice that for large 7 and not large ^7, there is an enormous "relativistic 
aberration": dt ~ dtsN/7^, and the observer sees a long CB story as a film in 
extremely fast motion. 

The energy of the photons radiated by a CB in its rest system, E^g, their energy 
in the direction 6 in the local SN system, Egj,^, and the photon energy E measured 
by a cosmologically distant observer, are related by: 

E^B = ^; El^ = {l + z)E, (3) 

with 6 as in Eq.(0). 



3 Reference values of various parameters 

To be explicit we must scale our results to given values of the parameters of the CB 
model. In this section we introduce the reference values that we adopt, which serve 
as benchmarks but imply no strong commitment to their particular choices. These 
values are listed in Table I, for quick reference. 
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Parameter 


Svnibol 


Value 


OIN oiidi o ilidoo 


Mo 

iVig 




SN-shell's radius 


Rs 


2.6 X 10^^ cm 


Outgoing Lorentz factor 


7out 


10^ 


CB's energy 


EcB 


10^^ erg 


Initial v^/c of expansion 


An 


1/(3 v^) 


Final /c of expansion 


Pout 


i/Vs 


Redsliift 


Z 


1 


CB's viewing angle 





10-3 



Table 1: List of the "reference" values of various parameters. In the text a barred 
parameter means its actual value divided by its reference value, so that, for instance, 
Ms = 1/2 means that the actual mass of the SN shell is taken to be SM©. Two 
parameters are not specific to our model (z and 6). 



Let "jet" stand for the ensemble of CBs emitted in one direction in a SN event. 
If a momentum imbalance between the opposite-direction jets is responsible for 
the large peculiar velocities of neutron stars, vns ~ 450 ± 90 km s~^ [^], the jet 
kinetic energy Ejet must be, as we shall assume for our GRB engine, larger than 
MnsVnsc ~ 10^^ erg, for Mns = 1-4 Mq []n|. We adopt a value of 10^^ ergs as the 
reference jet energy^. On average, GRBs have some five to ten significant pulses, so 
that the energy in a single CB may be 1/5 or 1/10 of Ejet- We adopt Ecb = 10^^ 
erg as our reference value. We denote with a bar the actual value of a parameter in 
the units of its reference value so that Ecb, for instance, means a given cannonball 
energy divided by 10^^ erg. 

Let 7in be the Lorentz factor of a cannonball as it is fired. We shall find 7in = 
(9(3 X 10^) to be a "typical" value (7in is not an "input" parameter). For this value 
and the reference CB energy, the CB's mass is very small by stellar standards, and 
comparable to an Earth mass: 



McB ~ 0.6 



3 X 10^ 

7in 



(4) 



The baryonic number of the CB is: 



E, 



CB 



2.2 X 10^^ Ecb 



nip c^ 7in 



3 X 10^ 

7in 



(5) 



The collision of a CB with a SNS is so violent — at ~ 1 TeV per nucleon — that 

^The jet-emitting process may be "up-down" symmetric to a good approximation, implying 
even bigger jet energies. In the accretion of matter by black holes in quasars 31 and micro- 



quasars ||2J] the efficiency for the conversion of gravitational binding energy into jet energy appears 
to be surprisingly large. If in the production of CBs the central compact object ingurgitates several 
solar masses, Ejet could be as large as ~ MqC^ ~ 1.8 x 10^^ erg. 
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there is no doubt that, as it exits the shell, the CBs' baryonic number resides in 
individual protons and neutrons. 

We have assumed that, in a SN explosion, some of the material outside the 
collapsing core is not expelled as a SNS, but falls back onto the compact object. For 
vanishing angular momentum, the free-fall time of a test-particle from a distance 
R onto an object of mass M is tfau = vr [R^/(8 G M)]^/^. For material falling from 
a typical star radius (R^, ~ 10^^ cm) on an object of mass M = 1.4 Mq, tfau — 1 
day. The fall-time is longer (except for material falling from the polar directions) 
if the specific angular momentum is considerably large, as it is in most stars. The 
fall-time is shorter for material not falling from as far as the star's radius. The 
estimate tfau — 1 day is therefore a very rough one. One day after core-collapse, the 
expelled SNS, traveling at a velocity vs ~ c/10 has moved to a distance: 

Rs = 2.6 X 10^^ cm (^) (1^) . (6) 

We adopt Rs = 2.6 x 10^"^ cm as our reference value. 

For the Lorentz factor of the CBs as they exit the SNS, we adopt the value 
7out = 10^, for the reasons discussed in the Introduction. Let /3inC be the expansion 
velocity of a CB, in its rest system, as it travels from the point of emission to the 
point at which it reaches the SNS, and let Pont c be the corresponding value after 
the CB exits the SNS, reheated by the collision. We expect these velocities to be 
comparable to the speed of sound in a relativistic plasma, c/a/3, as observed in the 
initial expansion of the CBs emitted by GRS 1915+105 [Q. As reference values, 
we adopt those of Table I. 



4 The collision of a CB with the SNS 



4.1 The shell's profile and transparency 



The density profile of the transparent outer layers of a SNS as a function of the 
distance x to the SN centre can be inferred from the photometry, spectroscopy and 
evolution of the SN emissions 
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The observations can be fit by a power law, x""^, 
with n ~ 4 to 8. Our results for neutrino fluxes are not sensitive to this density profile 
and our results for GRB 7-rays and for EGRs are only sensitive to the outer 
region where the SN shell becomes transparent. This implies that, for simplicity, we 
can adopt at all x > Rs the density profile observed in the shell's outer layers: 



p(x) =p(Rs)e(x-Rg 



X 



(7) 
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The SNS grammage still in front of a CB located at x is: 



roo 

Xs(x) = y p(y)dy = XsNS 



XsNS = 



4 7rRi 



(2.35 X 10^) 3 



Rs 

X 

Ms 

r; 



(8) 



For GRB photons in the MeV domain the attenuation length is similar, within 
a factor 2, in all elements from H to Fe, and it is close to the attenuation length 
in a hydrogenic plasma. In the CB model, at a fixed time, the energy spectrum in 



a GRB pulse is roughly thermal |22|. The radiation length in the obscuring shell. 



averaged over a black body spectrum of peak energy 1 MeV, is approximately: 



m 



>-GRB 



P ^ in „„v>.-2 



(TKN(lMeV) 



lOgcm-% (9) 



where ctkn is the Klein-Nishina cross section. For EGRs the attenuation length in 
hydrogen in the 100 MeV to 100 GeV range, dominated by e"*" e~ pair production, 
is O: 

XEGR^70grcm-2 . (10) 

The attenuation lengths of Eqs.(^ and (|10|) are all much smaller that the typical 
shell grammage of Eq. (P) . 

Equating Xs(x) and Xgrb and solving for x, one obtains the radial distance x^j^g 
at which the SNS becomes (one radiation length) transparent to GRB photons. For 
our reference parameters, some representative results are: 

Xgrb/Rs - (3.7, 6.2, 21) for n = (8, 6, 4). (11) 

The corresponding values for x^qj^, at a given n, are shorter: 

4gr/Rs - (2.3, 3.2, 6.9) for n = (8, 6, 4). (12) 

The GRB and EGR signals are emitted as the CB reaches the transparent outskirts 
of the SNS. The neutrino signal is emitted as soon as the CB starts colliding with 
the shell. We discuss in detail in Section 10 the time profiles and relative timing of 
these signals. 

4.2 Kinematics of a CB's collision with a SN shell 

The radius of the expanding CBs, as they reach the SNS, is: 

"3 X 10^ 



RcB ~ Rs — ^ 1.7 X 10^° cm A, 

7in 



7ii 



Rs , (13) 
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In its collision with the shell, a CB sweeps up a "target" mass 

Ms , (14) 



_ 

Mt ~ vtR^bXsns = Ms ^ 3.5 X 10"^ M^ Ar 

47in 



3 X 10^'^ 



7in 

where Xsns is the full column density of the shell, as in Eq.(^. 

Seen from the reference system in which the CB is at rest (and its shape, because 
of expansion, is roughly spherical) the constituents of the SNS impinge onto the CB 
with a Lorentz factor 7in. The average density of a CB with the reference radius of 
Eq.(p!3D and the reference mass of Eq.(^) is p ~ 1.8 x 10"'^ gr cm~'^. The nucleon- 
nucleon interaction length at that density is Acs = (Na c"™"*" p)~^ ~ 2.2 x 10^ cm, 
with Na Avogadro's number and o"™^ ~ 40 mb the proton-proton (or nucleon- 
nucleon) cross section at TeV beam energies. The CB's radius of Eq.(|13]) is much 
bigger than Xqb, implying that all nuclei in the region of the shell swept up by the 
CB interact. Approximately 1/3 of the energy in these collisions results in photons 
from tt" decay, which heat the CB to a temperature in the keV domain [0, |2z 



Seen from the reference system in which the SNS is at rest (or moving with a 
modestly relativistic velocity ~ c/10) a high-energy nucleon in the CB — suffering 
successive interactions in the dilute gas or plasma constituting the SNS — loses 
roughly 2/3 of its energy to vr^ production. The density of the shell is of order 
Ps = Ms/(4 7rRg) ~ 10^^^ gr/cm'^, for our typical parameters. At that density, the 
nucleon-nucleon interaction length is As ~ 5 x 10^^ cm, much less than the 0(Rs) 
shell's depth, so that the shell's material is, in this sense, "thick": it acts as a beam 
dump. The decay length of a charged pion of energy E^ is 5.6 x 10^ E^/(100 GeV) 
cm, much less than its interaction length, which is comparable to that of nucleons. 
Consequently, the beam dump is "thin" to vr decay and roughly 2/3 of a CB's nucleon 
energy is carried away by the neutrinos in vr — > p decays and in the subsequent p 
decays. 

The Lorentz factor of the CB after it has swept the SNS is simply the ratio of 
the total energy to the invariant mass of the outgoing object: 

Ecb/3 , - 

7out ^ , = , (15) 

'2 Mt c2 Ecb/3 + M^b 



where we have used Ecb ^ Mtc^, with Mt the target mass of Eq.([T^). Substituting 
for Mt and Mcb as functions of 7in and Pin, one obtains: 



/ 2 Eqb 

^ V 3AlMsc2 + 18EcB ^^^^ 



whose limiting values are: 



7out~^ (for6EcB»A'nMsc2) 



7out 



3 

7ii 



10 A, 



E, 



CB 



M. 



1 

2 

'2 A/r „2> 



(for6EcB<AnMsc^). (17) 
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For our reference parameters, Eq.(p!6|) implies that 7in ~ 10 7out- The very large 
"typical" values of 7in, ~ 3 x 10'^ or larger, as in Eqs.([T7|), imply that the fractional 
solid angle covered by a CB as it hits the SNS is tiny: An/(47ki) ~ 10~^ or smaller, 
for our reference jSi^. This presumably makes it unlikely for consecutive CBs to 
hit precisely the same spot in the SNS: CB-CB collisions and mergers may be the 
exception, rather than the rule, and the coUisional "histories" of successive CBs 
should be similar. 

Let (T^ ^ 6.5 X 10~^^ cm^ be the Thomson cross section, describing 7-e collisions 
at invariant masses comparable or smaller than the electron mass. The CB itself 
becomes transparent to the radiation it encloses when it reaches a radius R^g ~ 
[3McBCT^/(4 7rmp)]^ ~ 1.9 x 10^^ cm, for Mcb = 0.6 M^. If the CB in its rest 
system, after the collision with the SNS, is expanding at a transverse velocity /3out c, 
the distance away from the SN at which it becomes transparent is 7out R-cB/Z^out, or 
~ 3.3 X 10^^ cm for our reference parameters. By then, the CB is well out of the 
SNS and it has emitted from its expanding surface the radiation that constitutes 
the GRB signal [||]. 



4.3 Microscopic description of the collision 

The main point in outlining a microscopic picture of the collision of a CB and a SNS, 
as we shall see, is to conclude that the details of such a picture are immaterial to the 
estimate of the properties of GRBs and of their associated EGR- and high-energy 
u fluxes. But the discussion is important in that it sets the basis for how to make 
these estimates. 

Both the SNS and the CB are many pp interaction lengths long. The number of 
such lengths in the SNS is: 

NlTr^s = ^ - (5-6 X 10^) Ms [Rs]"^ , (18) 
TT rig 

As it enters the shell at a distance x = Rg from the SN centre, the number of pp 
interaction lengths in the CB is: 

'[AnRs]-'. (19) 

At a later point in the crossing of the SNS, e.g. at x = 2 Rg, when the CB is moving 
at 7 ~ 7out and expanding at a speed Pout, its number of interaction lengths is: 

N-* ~ (5.1 X 10^) EcB [7in]-' [Tout]' IPoM-' . (20) 

Thus, the number of pp interaction lengths in the CB is typically comparable to or 
bigger than the corresponding number in the SNS. 



N 



int 
CB 



Mcb in 
ttR^d 



A _TOT 



:i.6 X 10^) EcB 



7in 

.3 X 103 



10 



The simplest reference system in which to visuahze the colhsion is the centre-of- 
mass system (c.m.s.) of two slabs of nuclei, one belonging to the CB, the other to 
the SNS, both one interaction length long. Consider first the case in which the CB 
and the SNS are an equal number of nucleon-nucleon interaction lengths long. In 
the approximation of constant densities, the slab-slab c.m.s. coincides in this case 
with the overall c.m.s. of the CB and the material that it hits in the SNSQ. In this 
system both the CB and the shell are spatially contracted (relative to their respective 
rest systems) by the Lorentz factor ^Jy^jl at which their constituents are moving 
towards each other. After the nucleons have interacted once, their energy is degraded 
by an average factor f ~ 0.7, (the "leading particle" average-energy fraction observed 
in high-energy nuclear collisions). The nucleons of the leading slab, after the time 
required to interact a few times with a few "opposing" slabs, come to rest and are 
eventually turned back. Meanwhile fresh slabs are coming in and suffering the same 
fate as the first. An increasingly hot and dense pancake-shaped region is formed, 
containing the nucleons that have collided and the radiation initiated by 7's from 7r° 
decay. Because of its enclosed radiation pressure, this "pancake" eventually expands 
in its rest system at a speed of (9(c/-\/3). When all the slabs of the CB and of the 
SNS are consumed, the resulting object is the outgoing CB, at rest in this system. 

In the case where the number of interaction lengths in the CB and the SNS 
are different, a similar description applies in the system of reference in which the 
CB and shell densities are the same, up to the moment in which the object with 
the smaller number of interaction lengths is consumed. This object is typically the 
SNS, as Eqs.(|l8D to (pOD indicate. At that point, we have a hot and dense pancake- 
shaped object at rest, plus the fresh slabs from the CB's side of the collision that 
are impinging on the disk without having interacted yet. As these fresh slabs hit, 
they set the disk in motion. The final outgoing CB, now viewed from the local rest 
system of the parent SN, is moving with the Lorentz factor 7out of Eq . (p!6D . 

No doubt the previous description is oversimplified, for the violence of the colli- 
sion of the CB and the SNS presumably results in shocked and turbulent motions. 
Moreover, the freshly expelled SNS is probably not smooth, but also turbulent and 
inhomogeneous on small scales. Fortunately, a detailed description is not required 
for an estimate of the fluxes of z/'s and EGRs. 



■^It is easy to generalize the argument that follows to non-constant densities, by using a variable 
reference system in which the densities are instantaneously the same, but that requires an unjusti- 
fied amount of effort, the final results on the observable v fiuxes being the same. The generalization 
to unequal numbers of interaction lengths, we shall deal with explicitly. 
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5 The GRB 



In this section we briefly review the properties^ of GRBs, in the CB model that 
we need to establish comparisons between the GRB itself, and its accompanying 
EGR and v fluxes. 

In its rest frame, the front surface of a CB is bombarded by the nuclei of the 
SNS, which have an energy mpC^7 ~ 1 TeV per nucleon, roughly 1/3 of which is 
converted into 7-rays (from 7r° — > 77 decays) within a nucleon attenuation length: 



mr 



TOT 



42 gem 



(21) 



pp 



These high energy photons initiate electromagnetic cascades that, in turn, convert 
their energy to thermal energy within the CB. The radiation length of high energy 7's 
in hydrogenic plasma is Xegr, given by Eg . ([T0|) . The energy of the electromagnetic 



cascade ends up as heat, 
radiation length: 



The thermal photons, of energy E~ <C mgC , have a 



— ~ 2.0 g cm 



-2 



(22) 



The thermal energy contained in a CB's front-surface layer of "depth" Xt, con- 
tinually supphed by the SNS incident nucleons, will be radiated away. It is reason- 
able to expect that an equilibrium is established whereby, to a fair approximation, 
the quasi-thermal emission rate from the CB is in equilibrium with the fraction of 
energy deposited by the CB's collision with the SN shell within this one-radiation 
deep layer. A fraction Xx/Xp of the incoming protons interact in the radiating 
layer, and a fraction Xt/Xegr of the energy of the 7s from vr" production and de- 
cay is deposited in it. The total energy deposited by a single SNS shell nucleon is 
~ mpC^7/3. Equating the energy deposition per unit time to that re-emitted from 
the CB's surface as quasi-thermal radiation, we obtain an instantaneous tempera- 
ture: 



T(x) 





Xqi 


(n-l)c3[7(x)]21 


1/4 


X 


-n/4 




_Xp Xegr 






.-'^GRB. 





(23) 



where a is the Stefan-Boltzmann constant, 7(x) is a function evolving from 7i, 
to 7out5 XQf^g is as in Eqs.(0), and n is the SNS density index of Eq.(0). The 
temperature as the CB reaches the transparent region of the SNS, Ttp = T(xl 
is of order 0.15 keV, and is not very sensitive to the parameters of the model 
scaling roughly as: 

-n- 11 



tp 

GRB J) 

21, 



-tp 



oc 



Rr 



[7o 



11/2 



(24) 



since Xtp oc Rs and, in the outer regions of the SNS, 7(x) 



7out- 



We use the "surface model" of [E2| 
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The time-width of a single-CB GRB pulse is roughly characterized by a "trans- 
parency time" tcRB: the time elapsed between the moment the CB enters the SNS 
and the time it reaches its (one radiation length) transparent outer layer. In the 
observer's frame, this time is: 



l^GRB - 7 , (25) 

Tout c 



with Xqj^b as in Eg. (pJ]) . For our standard parameters and a typical 6 ~ 3/7out, 
tcRB is ~ 0.23, 0.45, 1.7 s, for SNS indices n = 8, 6, 4, respectively. 

The radius of the CB at the time t = tcRB is: 

Rc'b ~ RcB + (x^?RB - Rs) ^ (26) 

Tout 

with RcB as in Eq.(|T^). For times of C(tGRB), the radius of the CB increases 
approximately linearly with time. In the CB's rest system, and at a fixed time, 
the energy spectrum of the radiation emitted by the CB is an approximate black- 
body spectrum, corrected for absorption in the SNS, and emitted by a sphere whose 
surface grows as t^ and whose surface temperature decreases as [toRB/tj'''^''^^ as in 
Eg. (P^D . We have shown in |^ that this simple picture describes well the light 
curves and energy spectra of GRB pulses. The predicted GRB energy spectrum 
is the sum of thermal spectra with decreasing temperatures. Its high energy tail 
is exponential, with a characteristic temperature ~ Ttp. At high energies, this 
is an underestimate with respect to the observed GRB spectra, which decrease 
approximately as E~^. We attribute this discrepancy to the naivete of our thermal 
input spectrum: observations demonstrate that astrophysical plasmas subject to a 
flux of high energy particles — such as the CB in its rest system — radiate a "guasi- 
thermal" spectrum, corrected at high energies for such a power-law tail (clusters of 
galaxies are discussed in |]3^, galaxy groups in |^ and SN remnants in 



To estimate the total energy radiated by a CB's heated surface we must first 
compute the number Np of SNS nucleons that provide this energy, with the con- 
straint that the radiation they eventually produce be able to escape from the SNS. 
The naive estimate Np^^ ~ tt [Rcb]^ Xgrb/kIp turns out to be a very good approx- 
imation. Indeed, in terms of np(x), the nucleon number density in the shell, and 
XsNs = Xs(Rs), the shell's total grammage of Eg.(|), N^^^ is: 



GRB 



Xqrb 



/ dx TT [Rcb(x)] np(x) Exp -— ^ / np(x') dx' 

JRs L Aqrb "'X 

/■oo r, d r m 

/ dx [Rcb(x)]' —Exp --^ / np(x')dx' 

JRs ax L AoRB "'x 



vr [R^?b]' ^ (1 - Exp[-Xs/Xgrb]) ^ vr [R'^^f ^ , (27) 
mp mp 



13 



where we have approximated by a constant the radius of the CB as it travels through 
the outer few GRB absorption lengths. In the CB's rest frame, the total radiated 
energy is: 

prest _ XgrB X| TT [Rcb]^ 7out /p,3x 

%nlse ~ T? o • \^^) 

AeGR -^p -J 

As an example, for our standard parameters and n = 8, R^g ~ 4 x 10^^ cm and 
Ep'^uisc ~ 3 X 10^^ erg. The result scales roughly as R|/3oyt/7out and could be one 
or two orders of magnitude smaller for Rs and Pout somewhat below our reference 
values. 

An observer at rest, located at a known luminosity distance Dl(z) from the CB 
and viewing it at an angle 6 from its direction of motion, would measure a "total" 
(time- and energy-integrated) fluence per unit area: 



^, prest 



27out 



1 + 0' iLt 



(29) 



In a critical (17 = 1) Friedman universe the luminosity distance is given by: 



(l + z)c /-i dx 

where Hq is Bubble's parameter, Qu and Qa = Q — flu, respectively, are the matter 
and vacuum energy densities divided by the critical density (pc = 3HQ/87rG) and 
the radiation energy density has been neglected. In our explicit calculations we use 
Ho = 65 km/(s Mpc), = 0.3 and = 0.7, so that, for example, Dl(1) — 7.12 
Gpc ~ 2.20 X 10^^ cm. In Fig.(|) we show Dl(z) and [Dl(1)/Dl(z)]2 (the quantity 
to which we shall scale our results) for the quoted cosmology and, for comparison, 
for the case = 1, = 0, for which Dl(1) ~ 5.41 Gpc ~ 1.67 x 10^^ cm. 



6 The sources of high energy particles 
6.1 The origin of high energy z/'s 

The nuclei or nucleons in the incoming CB are comoving with it with a Lorentz factor 
7in. Nuclei in the outgoing CB have certainly been shattered into their constituent 
nucleons by the violence of the collision between the CB and the shell. They are 
comoving with the bulk Lorentz factor 7out- On average, a high energy nucleon 
colliding at a large centre-of-mass energy with a nucleon at rest — or moving along 
the same direction — exits the collision with a small transverse momentum and a 
fraction f ~ 0.7 of its original energy (for ultrarelativistic particles this statement is 
independent of longitudinal Lorentz boosts). The nucleons of the incoming CB must 
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have their Lorentz factor degraded from to 7out- On average, this takes a number 
i of high-energy coUisions satisfying f = 7out/7in5 that is i ~ 3, for ~ 3 7out5 or 
i ~ 6, for 7in ~ 10 7ouf 

Another way to reach a conclusion similar to the above is to view the interactions 
in the reference system introduced in Section 4.3. There we saw that the bulk of 
the CB's nucleons impinge, with a Lorentz factor ^/7J^/2, on a "pancake", which 
is at rest, or moving towards them at a mildly relativistic velocity ~ c/v^- For 
7in = 3 X 10^ (10^) the incoming energy of the CB's nucleons is E ~ 26 GeV (~ 47 
GeV). Consider the number of interactions i necessary to bring these nucleons down 
to an energy (~ 5 GeV) below which the multiplicity of pion production on a 
stationary target is no longer roughly constant (up to logarithmic corrections), but 
is suppressed by threshold effects. The argument of the previous paragraph now 
yields i ~ 4 (6) for 7in = 3 x 10^ (10^). One may be concerned with the fact that, 
if the CB contains more interaction lengths than the SNS's target funnel, the last 
of the CB's protons to interact with the pancake (which is by then moving in the 
CB's direction) may only suffer collisions at a centre-of-mass energy insufficient to 
produce pions. To dissipate this concern, consider the very last nucleon of the CB 
to suffer a collision with the ensemble of the CB plus the swept-up mass, and go 
back to the system in which the SN is at rest. In that system, the incoming nucleon 
and its target have Lorentz factors 7in and 7out, respectively, which brings us back 
to the argument in the previous paragraph. 

Another concern is that, if nucleons suffer only a few interactions with a large 
centre-of-mass energy, the earlier estimate that 2/3 of the CB's energy is lost to 
neutrinos may be grossly incorrect. But after i interactions the fraction of the 
original energy of the CB's nucleons that has been invested in pion-pro duct ion is 
a ~ 1 — f, two thirds of which ends up in neutrinos. For i = 3 (6), a ~ 66% 
(~ 88%), so that Eq.(|TBp is a fair approximation. 

Given the previous arguments, we shall estimate the neutrino flux as that pro- 
duced by a nucleon beam, containing as many nucleons as the incoming CB, moving 
with a Lorentz factor 7in, and interacting thrice on a nucleon or nuclear target. Be- 
cause of "Feynman scaling" , as we shall see, it does not matter whether the target 
nucleons are at rest or receding with a Lorentz factor 7out- It follows from the above 
discussion that this estimate of the neutrino flux is an underestimate. However, in 
practice the incurred error is small, because the charged pions produced in i > 4 
interactions have a relatively low energy. Since the detection efficiency of their decay 
neutrinos is weighted, as we shall see, by two powers of energy, the low-energy tail 
of the neutrino spectrum is immaterial. 
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6.2 The origin of EGRs 



The SNS is only transparent to 7-rays in its outer layer, some Xegr = 70 g cm~^ 
deep, Eq.(|l^). That figure also corresponds to roughly two high-energy nucleon- 
nucleon interaction lengths, see Eg. (piD . The shell is many interaction lengths thick, 
so that by the time the CB reaches the shell's 7-ray transparent outer layer, it is 
already moving with 7 ~ 7out. At that point the CB has reached a radius: 



R, 



CB 



(^EGR 



but 



7out 



(31) 



For Pout = l/v3, 7out = 10 and the largest of the x^qj^ in Eg. ([1^) , Rcb ~ 8.8 x 
10^^ cm. With this radius and the reference mass of Eg.(^) the CB's grammage is 
Mcb/I^tRcb) ~ 1500 g cm~^, which is larger than that of the transparent outskirts 
of the shell, Xegr- We shall therefore compute the EGR fiux as that originating 
from the vr^'s made by the front of the CB as it interacts with all the nucleons in 
the shell's transparent outer layerQ. This total number, computed in the same way 
as N^RB in Eg.(^, is: 

Xegr ^32) 



N 



EGR 



vr R; 



whose numerical value is: 
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EGR 
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(33) 



where the values of x^qj^ for various shell density indices n are those of Eg. ([T^) and 
we have not made explicit the weak dependence on (Rg/Mg) to the power l/(n — 1). 

The total energy of the EGR pulse generated by a CB, as seen by a local observer 
at rest in the SN system, is: 



E 



1 



EGR 



TT Rcb Xegr 7out • 



(34) 



To study the spectrum and angular coUimation of the EGR fiux, as seen by a 
cosmologically distant observer, we must recall the details of pion production and 
decay. 



6.3 Pion production in nucleon— nucleon collisions 

The CB's baryon number, as it crosses the SNS, resides in protons or nuclei that are 
been broken into their constituents by relativistic collisions. The same is the case 

^We neglect the reinteractions with the shell's nucleons that have already been struck one or 
more times and set into forward motion, thereby slightly underestimating the EGR flux. 
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for the funnel in the SNS that is swept up by the CB. At TeV energies, the nucleon- 
nucleon, nucleon-nucleus or nucleus-nucleus processes have different cross sections, 
but the properties of the produced pions, per colliding nucleon-nucleon pair, are 
very similar, and not significantly different for protons or neutrons. Consequently, 
we can use for our considerations the empirical information on pion (and kaon) 
production in proton-proton collisions. 



Bailly et al. |^ reported on a study of inclusive charged pion production in the 
collisions of protons of energy 360 GeV on a hydrogen target (c.m.s. energy ^/s ~ 26 
GeV), function of "Feynman x" : 

x.^. (35) 

The result, roughly the same within errors for Tr"*" and vr" is: 

F.W-^/dp;^^0,2^(l-x)3.. (36) 

where p.^, is the pion transverse momentum. Given the approximate isospin inde- 
pendence of the interactions, the above result should also apply to the inclusive vr" 
production. The hypothesis of Feynman scaling, satisfied up to small logarithmic 
corrections, is that Eg. (P5|) is independent of energy. 

The X and p^ dependences are observed to factorize to a good approximation, 
and the dcr'^/dp^ distribution is roughly exponential, with average 

Pl ~ 320 MeV. (37) 

The most precise data on transverse-momentum distributions at the relevant c.m.s. 
energies were collected for 7r° production in pp interactions at the ISR collider PB[; 



the measured single-photon yield resulted in p^ ~ 160 MeV, whence the result we 
adopt for neutral or charged pions (p^ ~ 2 p^ ) . The double differential cross section 
for inclusive vr production is therefore of the form: 



da^ 



aJpOT dp2 dx 



F^(x) G^(p^ 



G.(pJ^ -l^e'PT/Px, (38) 

Z Prp 

with F^(x) as in Eg . (|36|) and p.^, as in Eg. (|37D . For ultrarelativistic reaction products, 
Eg.(p8D is invariant under longitudinal Lorentz boosts, so that, with x ~ Ejf'^/Ejf'^, 
it can be used for proton interactions on a stationary target. We shall use it for 
incident Ep up to a few TeV. 
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7 The flux of EGRs 



To compute the spectrum of outgoing photons per nucleon-nucleon collision, we 
must convolute Eq. (|38D with the distribution of tt^ 77 decay. For ultrarelativistic 
pions, the distribution of fractional photon energies (w = E^/Ett) is fiat and limited 
by < w < 1. The 7 distribution in y = E^/Ep, produced by the decay of pions 
distributed as in Eg. (PUD , is: 



F,(y) = 2 F.(x) ^ dw 5 fw - -] , (39) 



where the prefactor is for the two 7's per ir^ decay. To a few per cent accuracy, the 
result of the convolution can be fitted^ by: 

F,(y) ^ A,ie-^-y (40) 
~ 1.1 , b^ ~ 8 . (41) 



The 7-production double differential cross-section in y and p'*' is of the same form 
as Eq . (|38|) , with a longitudinal factor F^(y) and a transverse factor with p''' ~ 160 
MeV. Let E^ be the energy of a photon as it reaches the Earth, cosmologically red- 
shifted by a factor 1 + z, and let Ep ^ mp c^ 7out be the energy of the CB's nucleons, 
in the local rest system of their SN progenitor, as they reach the outer part of 
the SNS. For the small angles 6 at which the 7-rays are forward-collimated by the 
relativistic motion of the parent tt^'s, the photon-number distribution in x^ = E^/Ep 
and cos 6', per single nucleon-nucleon collision, is: 



dn^ 
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(42) 



Let dn^/dQ be the total (time-integrated) number flux of EGR photons per unit 
solid angle about the direction 6 (relative to the CBs' direction of motion) at which 
they are viewed from Earth. The photon number distribution per incident CB is: 



dn- 



dx^ dQ 



N 



EGR 



B. 



f7(z,7out,6',x^) 



~ x^ e 



(43) 



^An exponential fit is inadequate close to the limit y = 1, but for y > 1/2 the flux is negligible. 
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with Np*^^ given by Eq.(^). Since a typical GRB has an average of ncB = 5 to 
10 significant pulses, the total fiux of EGRs in coincidence with a GRB may be an 
order of magnitude above that of Eq. f^Sf ). In Fig.(^) we show f^ as a function of 
at various 6; for z = 1 and 7out = 10^. The average fractional EGR energy in 
the spectrum of Eq.(|^) is = 2/c^, corresponding, at z = 1 and for 



average energies E-y ~ 120 GeV for ^ = 0, 



E^ ~ 70 GeV for 6 = l/7out, and E^ 



10^ to 
40 

GeV for 9 = 3/7out, a more probable angle of detection [Q. Except at the highest 
of these energies and/or at redshifts well above unity, the absorption of 7-rays on 
the infrared background — for which we have not explicitly corrected Eq.(^) — is 
negligible. 

Roughly characterize the efficiency of a 7-ray detector as a step function 0(E'^ — 
The total fiux above threshold, per incident CB, is then: 
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(44) 



where the scaling properties of Np*^^ are those of Eq. (|33|) . In Fig.(^) and @d) 
we show Gy as a function of xj^jj^ at various fixed 9, and vice versa; for z = 1 and 
7out = 10^. The very large fiux dn!!l^[0, 0]/df2 of Eq.(^4|) is seen to be significantly 
reduced as soon as 9 and/or x^j^ depart from zero: the EGR fiux is not as gigantic 
as it appears to be at first sight. 



7.1 EGR versus GRB total energies 

An observer at a cosmological distance from a GRB source, if equipped with a 
detector of sufficient angular coverage, would measure a total energy per CB, in the 
MeV-range photons of a GRB pulse, of E^j^g = 7out Epufse/ (1 + with Ep'^^^gg given 
by Eq . (pSD . In the multi-GeV range of the EGRs, the measurement would result in 
E|gr = Eegr/(1 + z), with E^^^^ given by Eq.(|4|). The ratio of EGR to GRB total 
energies is: 
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(45) 
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The factor l/7out may look surprising at first: the total GRB energy is oc 7^, while 
the EGR energy is oc 7. In both cases, one factor of 7 is associated with individual 
nucleon energies. But GRB photons "benefit" from the energy associated with the 
bulk motion of the CB, which acts as a relativistic mirror (target SNS particles at 
rest, if elastically and coherently back-scattered by the much heavier CB, would 
recoil with a Lorentz factor 7^, while particles produced in individual pp collisions 
would carry an energy scaling as 7). 

Using the results of Eqs.(|lT|) and (|I2|) for the transparency distances, we obtain: 

^ (0.71, 0.54, 0.26) for n = (8, 6, 4), (46) 

that is, the EGRs carry almost as much energy as the GRB. Since the individual 
EGR photons have energies of C(IOO) GeV, as opposed to 0{l/2) MeV for GRB 
photons, the number of EGR photons is five or six orders of magnitude below that 
in the associated GRB. 



8 The flux of high energy neutrinos 

The calculation of the flux produced in the collision of a CB with the SNS is 
analogous to the calculation of the photon flux. The z/^ flux gives rise to a signal of 
about 1/3 the size of that of the flux (we neglect it, since we find it preferable 
to establish a lower limit to the observational prospects). The i/^'s are made in 
the chain reactions pp — tt + Tr"*" /W"^ t'^; and tt^ i^n, followed by 

/i" e~ Vf.. We have also estimated the contribution of K production and decay, 
which turns out to be negligible. 

8.1 Pion and muon decay 

We have shown in section 6.1 that, in order to estimate the neutrino flux from 
the CB-SNS collision, it is adequate to consider an average of i ~ 3 interactions 
of the incoming nucleons. The leading outgoing particle in these interactions has 
an average transverse momentum significantly smaller than that of the produced 
mesons: it can be neglected. The simplest way to compute v fluxes is to work out first 
the pion distributions made in three successive interactions of the incoming nucleons 
and then convolute the result with the pion decay distributions. The pions made in 
the first interaction have the distribution of Eq.(^6|). To compute the distribution of 
those made in the two successive interactions, we must make convolutions, analogous 
to that in Eq . (|39D , with the distribution of the exiting leading- nucleon longitudinal- 
momentum distribution. For the latter, we use the measurements of Ref. [^. The 
result, of summing the pion longitudinal distributions as functions of x = E^r/Ep 
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(with Ep the original incoming nucleon's energy) can be simply parametrized, to 
~ 15% accuracy, as: 

F:'(x)^[1 + 2.2(1-x)^]F.(x), (47) 
with F,r(x) given by Eq.(p6|). 

In the decay in flight tt^ fi^ of pions with E,r ^ hItt, the distributions of 
fractional neutrino and muon energies (x;^ = Ei,/E,r and x^ = E^/E,r) are flat and 
limited by < x,^ < x^ax = 1 — m^/m^ and 1 — x^ax < < 1- The z/^ distribution 
in y = E,^/Ep, produced by the decay of pions distributed as in Eg . (p6|) , is given by: 

F.(y)=/"F:'(x)*r"^^(x„-^) ; (48) 

Jo X Jo Xinax V X/ 

the muon distribution in the decay vr" — > /i~ is analogous, with the proper change 
of integration limits. The distribution in z = E,^/E^ in the decay of left-handed 
muons is 3z^ and, upon neglect of m^/m^, it extends from to 1. We do not 
write here explicitly the double convolution, analogous to Eg . (|48|) , involved in the 
calculation of the y-distribution in vr — > — ^ z/^ decay. 

The mean Xj, = E^^/Ett in tt —>■ fiu decay is x^, = Xniax/2 — 0.19, while for the 
muon ~ 0.81. The mean z = E,^/E^ in vr ^ z/^ ... decay is z = 3/4, so that the 
mean x^, in the tt — > yU — > z/^ decay chain is ~ ^ ~ 0.6. The available rest energy 
in TT — > /xz/ or /i — i> euu decay is small relative to the mean transverse momentum of 
the parent pion. This implies that the mean z/^ transverse momentum in the vr — /xz/ 
chain is x^, px — 60 MeV, where we have used Eg.(^7|). The corresponding result for 
the TT — i> /X ^ z/^ chain is px — 190 MeV. We shall see that the muon detection 
sensitivity on Earth is weighted by two powers of energy (one for the cross section, 
one for the muon range), so that the vr ^ /i ^ z/^ process, which produces a harder 
z/^ beam, is harder. Rather than giving results for a two-component distribution 
(vr fiu and /i euu) we shall use a common transverse momentum: 

p"^ = 190 MeV (49) 

for the overall z/^ beam. This results in a small underestimate of the flux at a fixed 
angle. 

The relative contribution of kaons to the z^^ flux is suppressed with respect to 
that of pions for three reasons. The K/vr relative multiplicity is ~ 1/5; the K — > /xz/ 
branching ratio is ~ 63%; and the available energy in the decay is not negligible 
in comparison with p^. All in all, the K-decay contribution to the z/^ flux at fixed 
angle is at the few per cent level: we neglect it altogether. 

The result of all this analysis is a longitudinal distribution in y = Ej,/Ep (with 
Ep the incoming nucleon's energy) that can, to a few per cent accuracy, be fitted 
by: 

F,(y) ~ A.-e-^-y (50) 
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A. 



12 . 



(51) 



The z/^-production double differential cross-section in y and p';^ — describing the 
neutrino flux generated in the beam dump per incident proton — is of the same form 
as Eg. (PSP , with a longitudinal factor Fy(y) and a transverse factor with p';^ ~ 190 
MeV. 



Let Ej, be the cosmologically redshifted energy of a neutrino as it reaches the 
Earth, and let Ep ~ mp c^ 7in be the energy of the CB's nucleons, in the local rest 
system of their SN progenitor, as they enter the SNS. In analogy with Eq.(|^) the z/^- 
number distribution in x,y = E^/Ep and cos 6, per single nucleon-nucleon collision, 
is: 
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(52) 



Let dn^/dQ be the time-integrated number of neutrinos per unit solid angle 
about the direction 6 (relative to the CBs' direction of motion) at which they are 
viewed from Earth. In analogy with Eq. (|43|) , the neutrino number distribution, per 
incident CB, is: 
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(53) 



with Nb the total baryon number of the CB, given by Eq.(^. For a GRB with ucb 
signiflcant pulses, the total number of neutrinos is ucb times larger than that of 
Eq.(PD. 



In Fig.(PD) we show f^ as a function of x^ at various 9] for z = 1 and 7in = 
10"^. The average fractional v energy in the spectrum of Eq.(^) is x^ = 2/c^, 
corresponding, for the chosen z and 7in, to average energies E^ ~ 712 GeV for = 0, 
E^ ~ 315 GeV for 9 = 1/10^, and K ~ 150 GeV for 9 = 3/10^. 

Neutrino oscillations may reduce the flux of z/^s of Eq.(p3D by as much as a factor 
of 2 (if they are maximal) or even 3 (if they are "bimaximal" ) . 



8.2 Muon production on Earth 



Muon neutrinos produced by a GRB can be detected by large-area or large-volume 
detectors, in temporal and directional coincidence with a GRB 7-ray signal. The 
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detection technique typically involves the "upward-going" muons, for which there is 
no "atmospheric" cosmic-ray background. 

A flux of neutrinos traversing rock or ice interacts with target nuclei N, producing 
muons in the process z/^+N — > fi+... In the energy range of interest here, the inclusive 
muon cross-section per target nucleon is: 



dcr(E^, E^) 
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(54) 



The produced muons lose energy and "range-out" in matter before they decay. 
At the energies of interest here, the muon energy loss per unit distance x (in a 
material of average atomic number and mass Z and A) can be approximated by: 



_ ^ = R(E)^-^Ro(l + BE) 
dx pw 

where p is the material's density, pw is 1 g cm~^, and 
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(56) 



In ice or a typical rock material 2 Z/A ~ 1, and Z is small enough for the neglect of 
the B term in Eg . (|55|) , at the muon energies we shall encounter (E^ ^ 1 TeV), to 
be a good approximation. 

At a given position x in a target material, an (approximately x- independent) 
flux per unit area dNi,(E^) / dEj, dA gives rise to a p flux dN^(E^, Ej,, x) / dEj^ dE,^ dA 
satisfying the equation: 



d_ 
dx 



dN, 



dE„ dE„ dA 



pNa 



da dN^ 



dE^ dE^dA 



dE d 



dx dE,. 



dN. 



dE„ dE„ dA 



(57) 



For a target thickness much larger than the muon range, an equilibrium between the 
produced and slowed-down muons is reached, whereby the muon flux is independent 
of position and the l.h.s. of Eq. (^7D vanishes. Inserting Eqs. (|53|) and ( |5^ ) into Eq.(^7D 
and integrating, we obtain ||40[| : 



dN,, 



dE^ dA 



dN 

K dE. (E. - E, 

E„ dE^ dA 



K 



n.0 Jii. 



2.26 X 10-^^ GeV 



(5^ 



Define x^ = E^/Ep: the ratio of the energy of a muon produced on Earth to the 
energy Ep = mpC^7in of the CB's nucleons, as they enter the SNS. Substitute the 
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neutrino flux of Eq. ( pSf ) into Eq. 
a muon flux per incident CB: 



dn„ 



dx^ dQ 



KE 



and integrate over neutrino energies to obtain 
dn,. 



P A„ dx, dl] ^""^ 



KE^|^(l + ^rf.(z,7in,^,x.) 

n3 ^ 



(59) 



with and as in Eq.(p^ and Nb the total baryon number of the CB, Eq.(^. In 
Fig. (0c) we show f^ as a function of x^ at various 6, for z = 1 and 7in = 10^. 

Very roughly characterize the efficiency of an experiment as a step function 
jumping from zero to unity at E^ = E[^jj^. The observable number of muons per CB 
and per unit area, obtained by integration of Eq.(^), then is: 



dnT[x:;,„,^] dnT[0,0] 



G, 



dn 

;i+z)b. 



GM(2'7in,6',x; 



ram J 



1 + 



g '^•^ ^min 



dn^[0,0] 
dO 



E^. 

mm 

nip7in 
3.2 X 102_ 
km^ ""^ 



7ir 



104 



Dl(1) 
Dl(z) 



(60) 



In Figs.(^,b) we show as a function of x(^jj^ at various fixed 6, and vice versa; 
for z = 1 and 7in = 10^. The relatively large flux dn^[0,0]/dn of Eq . (|60D is seen 
to be very significantly reduced as 6 and/or x(^j^ depart from zero. Once again, for 
a GRB with ucb significant pulses, the total number of muons is ucb times larger 
than that of Eq . (|60|) , and neutrino oscillations may reduce the flux by a factor 2 
or 3. 



9 Angular apertures and observational prospects 

Barring the case of GRB 980425 — whose exceptional properties and their interpre- 
tation within the CB model are discussed in |^ — the equivalent spherical energies 
of the GRBs with measured redshifts range between ~ 2 x 10^^ erg (GRB 990123) 
and ~ 2 X 10^^ erg (GRB 970228). The dependence of the GRB flux on the angle 
9 subtended by the CB's velocity vector and the hne of sight is given by Eq.(^): 
df/dfl oc (1 + 6'^7o^J^'^. This 9 dependence is the steepest parameter dependence 
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of the CB model, see Fig.l of ||22[. It is therefore reasonable to attribute the range 
of observed equivalent spherical energies to the 9 dependence, as if GRBs were oth- 
erwise approximately standard candles. The observed three orders of magnitude 
spread in equivalent energy then corresponds, according to Eq.(^), to a spread of 
viewing angles between 6' ~ and 6 ~ 3/7out- The cutoff at the upper angle reflects 
the sensitivity of past and current observations. 

The energies of the individual GRB 7-rays and the GRB fluences indicate CB 
Lorentz factors 7out ~ 10'^. So does an approximately 1:1 SN/GRB association. 
(The geometrical fraction of currently observable GRBs, for ^ < 3/7, is 7r^^/(47r) 
9/(47^). For 7 = 10'^ this fraction precisely reconciles the SN II, lb, Ic rate in the 
observable universe: 12 ± 5 s~^ [41] with the corresponding GRB rate of ~ 10^ per 
year.) 

The fluxes of i/-induced muons and of 7's of GRB and EGR energies have different 
9 dependences and the circumstance that GRBs are currently observed at angles up 
to ^ ~ 3/7out ~ 3 X 10~^ plays an obvious role in discussing the search for EGR 
and V signals in spatial coincidence with GRBs. The discussion is summarized in 
Figs.(^), where we compare the angular apertures of the three fluxes. The absolute 
and relative normalizations in these figures are arbitrary, so that the GRB results, 
based on Eq. , depend only on 7out, chosen to be 10^. The EGR results, based on 
the second of Eqs. (^41) , depend also on z (chosen at z = 1) and on Ej^j^, taken here 
to be 50 GeV. The v results, also for z = 1, are based on the second of Eqs.(pll|); 
they are for 7in =10"^ and E^^j^ = 50 GeV. 

According to Figs.(^, the EGR beam, up to very large 9, has a broader tail 
than the GRB beam. In practice that means that a detector with the sensitivity 
to observe the EGR flux of Eq.(|^) should find a signal in temporal and angular 
coincidence with a large fraction of detected GRBs. The i/^-induced beam is 
about an order of magnitude narrower than the GRB beam in angle, two orders of 
magnitude in solid angle. Consequently, a detector with a sensitivity close to that 
necessary to observe the /i flux of Eq.(^) would see coincidences with only about 
one in a hundred intense GRB events, that is ~ 1% of GRBs in the upper decade 
of observable fluences, for which 9 ~ l/7out- 

To ascertain the observational prospects for EGRs and i/'s, one would have to 
convolute our predicted fluxes with the sensitivities of the many large-area or large- 
volume V and EGR "telescopes" currently planned, deployed or under construction. 
We do not have sufficiently detailed information to do so, but a coarse look at their 
potential indicates that testing the CB model will neither be trivial, nor out of the 
question. The small area of past detectors with a capability to see EGRs, such as 
EGRET, would preclude the observation of the flux of Eq.( 
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10 Timing considerations 



In the cannonball model, each CB crossing the SNS generates an individual 7- 
ray pulse in a GRB light curve. The complementary statement need not be true: 
not every observed pulse necessarily corresponds to a single CB, since the 7 rays 
generated by sufficiently close CBs may overlap. This can be seen in the two top 
entries in Fig.(|^), which show the lightcurves of the same ensemble of CBs crossing 
two SNSs, which differ only in their density-profile index; in the case of the more 
extensive SNS (n = 4) the various CBs blend into a single pulse. 

Each CB should generate three distinct pulses: a GRB pulse, a v pulse and an 
EGR pulse. The v and EGR pulses are narrower in time than the GRB pulse and 
they preceed it. Observed with neutrinos or EGRs, then, a burst has the same pulse 
structure as the GRB, but the pulses are shorter and are precursors of the GRB 
pulses. We proceed to estimate the magnitude of these effects, illustrated in Fig.(0). 

For a given density distribution of the SNS, such as that in Eq.(^, it is possible, 
though laborious, to explicitly compute the expected time profile of the neutrino 
signal. This proffie is sensitive to the shape of p(x) at all x, including the inner part 
of the shell, for which no empirical data are available. Consequently, we shall only 
give here approximate results for the width in time of the v signal, and for its timing 
relative to the onset of a GRB pulse. 

Let be the distance from the SN centre at which the shell's grammage, as in 
Eq.(|D, is half of the total SNS grammage, that is = Rg 2^/*^^^"^). The temporal 
half-width of the neutrino signal, tjy, is roughly the time it takes the CB to reach 
this pointQ. As measured by the observer, this time is given by the same expression 
as Eq.(^) with the substitution of Xqj^b by x^. The ratio of durations of a single 
pulse in neutrinos and in GRB 7-rays is: 

ti, x,^ — Rs , , 

; 15 (61) 

tCRB X^jj^B - Rg 

For SNS density indices n = 8, 6 and 4, this ratio is 0.038, 0.029 and 0.013, respec- 
tively: the duration of a pulse is a few per cent of that of an individual GRB 
pulse. Neutrinos are emitted from the moment the CB hits the SN shell, while GRB 
7 rays can only be seen if emitted from the transparent SNS outer layer. The time 
difference between the onset of the corresponding pulses is ~ tcRB: a z/ pulse should 
precede its corresponding GRB pulse by approximately the width of the GRB pulse. 



The discussion of the EGR pulse follows analogous lines. The ratio of EGR and 
GRB pulse widths is: 

tsGR Xegr - Rs 



tcRB X^J'rb - Rs 
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^The time needed for the "last proton" of the CB to catch up and interact with the rest of the 
colhding CB is shorter than t^ by a factor ^ /3in (7out/7in)^ R-s/(x;/ — Rs)- 
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with x^GR given by Eq.(|l^). For SNS density indices n = 8, 6 and 4, this ratio is 
0.48, 0.42 and 0.29, respectively: the duration of an EGR pulse is shorter than that 
of the corresponding GRB pulse by a factor 2 or 3. The time difference between the 
onset of the corresponding pulses is a fraction 1 — tECR/tcRB of the duration of the 
GRB pulse: an EGR pulse should preceed its corresponding GRB pulse by 50 to 
70% of the width of the GRB pulse. 

The light curve of an EGR pulse is proportional to the SNS shell density, cor- 
rected for absorption. For the density profile of Eq.(0) and the corresponding gram- 
mage of Eq. (PD : 



dN 



EGR 



dt 

x[t] 



oc P3(x[t]) Exp 



Xs(x[t]) 



X 



EGR 



1+Z 



Ct 



(63) 



The considerations of this section are visualized in Fig. (0) , where we have drawn 
the light curves of a single GRB in GRB 7-rays, in EGRs and in neutrinos. The 
timing sequence of the pulses is put in by hand and their normalizations correspond 
to (random) values of 7out close to 10^, see Eq. (p8|) . The two columns of the figure 
correspond to n = 8 and n = 4. Notice how the EGR pulses precede the GRB pulses 
and are narrower: the EGR has a better time "resolution". For neutrinos, this is 
even more so. 



11 Conclusions 

In the CB model of GRBs, illustrated in Fig.(|I]), cannonballs heated by a collision 
with intervening material produce GRBs by thermal emission, and their electron 
constituency generates GRB afterglows by bremsstrahlung, synchrotron radiation 
and inverse Compton up-scattering of these photons and the cosmic background 
radiation. The material CBs hit is an excellent "beam-dump", so that nucleon- 
nucleon collisions generate a very intense and coUimated flux of neutrinos. Because 
of absorption, the emission of energetic 7-rays via vr" production and decay is much 
less efficient, but by no means negligible. 

The u flux has a total energy of the order of 10^^ erg (roughly 1/3 of the total 
energy in a jet of CBs, augmented by the ratio 7out/7in5 and reduced by the redshift 
factor). But individual neutrinos have energies of only a few hundred GeV, as 
illustrated in Fig.(|^), and their enormous flux will be hard to detect, even though it 
is collimated within an angle ~ 10"'^. The detection in coincidence with GRBs will 
be further hampered by the fact that the GRB angular distribution is broader, as 
shown in Fig.(|D. 
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The EGR flux carries roughly as much energy as the GRB, that is Epyfgg7out ~ 
]^q48 gj^g pgj, puise^ with E^p-^is^ as in Eg. (p9|) . The EGR beam, as shown in Fig.(||), is 
somewhat broader than the GRB beam, so that the search for coincidences should be 
fruitful. The typical energies of EGRs, as illustrated in Fig.(||), are of tens of GeVs, 
and the relatively high threshold energies of current large-area detectors should be 
a limiting issue, as in the case of neutrinos. 

The pulses of the GRB 7-rays should be slightly preceded by narrower pulses 
of EGRs and by much narrower pulses of i/'s, as illustrated in Fig.(j^). The CB 
model, as we have seen, predicts very specific properties and relations between the 
GRB, EGR and u spectra and light curves. In this respect, as in many others, the 
Cannonball Model is exceptionally falsifiable. 
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Figure 1: The CB model in a SN environment, not shown to scale. Relativistic CBs 
are emitted by a compact object accreting matter from a disk and/or torus. They 
hit a SN shell generating z/'s, quasi-thermal radiation (the GRB) and 7-rays from 
vr*^ decay (the EGRs). The latter two exit only from the transparent outer layers of 
the SN shell. 
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Figure 2: Luminosity distances and ratios thereof, as functions of redshift, for two 
0=1 Friedman universes, with two choices of matter and vacuum densities. 
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Figure 3: EGR, neutrino and muon fluxes, at various fixed observation angles 9, 
as functions of the fractional momentum of the observed particle, at redshift unity. 
The functions f^(z, 7out, 6*, x^) of Eq.(m), for 7out = 10^ and f^,^(z, 7out, 6*, x^,/.) of 
Eqs.(|53|, |5^) , both for = 10^, are depicted. 
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Figure 6: Comparisons of angular distributions of GRB photons, EGR photons and 
i/-produced muons in water or ice. In the upper graph, the normahzations of the 
three curves are arbitrary. In the lower one, they are all normalized to unity at 
^ = 0. 
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Figure 7: A "synthetic" 7-ray burst consisting of five CBs with 7out within a factor 
of 2 of 7out = 10^, with other parameters at their reference values. The CBs are fired 
at random times in a 2.5 s intervaL The two columns are for SNS density indices 
n = 8 and 4. Top: the event seen in the 30 keV to 1 MeV GRB domain. Middle: 
seen in EGRs from tt^ decay. Bottom: the neutrino signal. 
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